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Abstract 

Analyzing the production of weakly relativistic plasma produced by mi- 
crowave fields with circular polarization the electron distribution function is 
obtained to be non-equilibrium and anisotropic. Furthermore, it is shown that 
produced plasma is accelerated on the direction of propagation of microwave 
electric fields. The electron velocity on this direction strongly depends on 
electron origination phase, electric field phase, and amplitude of microwave 
electric field. Making use of the dielectric tensor obtained for this plasma, it 
is shown that the wiebel instability develops due to the anisotrpic property of 
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distribution function. The dispersion equation is obtained for this instability 
and the growth rate of it, is calculated. 
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I. INTRODUCTION 



A microwave produced gas discharge is a rather complicated phenomenon ex- 
hibiting a variety of features. Numerous theoretical and experimental studies have 
been devoted to this phenomenon. The interaction between intense (MW) fields 
and neutral gas open new possibilities for studying the fundamental property of this 
produced plasma. In a plasma produced in the interaction of the MW field with a 
gas, the electron distribution is non-equilibrium and may give rise to various plasma 
instability. 1 The main features of this interaction is the gas ionization mechanism. 
Furthermore, the key role in the plasma processes is played by the kinetic effects 
associated with the specific features of the electron distribution. Due to the short 
interaction times, these features are governed completely by the pulse parameter of 
microwave field. In addition the interaction between intense MW fields with an inho- 
mogeneous plasma results in variety of phenomena such as frequency up shift during 
the propagation of a pulse microwave field through a plasma. 2 In such a strong wave 
field the electron oscillatory energy e e is much higher than the ionization energy Ii on iz 
of gas atoms 

e 2 El 

where u is the radiation frequency, m is electron mass, and E is the electric field 
strength. When field amplitude is comparable to the atomic field strength E a ps 
5.1 x 10 9 VcvrT x the tunneling ionization becomes an important mechanism for di- 
rect ionization of the gas atoms. This effect was completely studied in the different 
papers. 3,4 ' 5 On the other hand relativistic effects come in to play when kinetic energy 
of electron oscillation in an electromagnetic field is comparable with the electron rest 
mass energy. Moreover, the microwave fields generated by the present day pulsed 
duration sources are weaker than the atomic field and are capable of manifesting 



weakly relativistic behavior of electron produced. On the oder side, the aforemen- 
tioned effect can be easily manifested by intense laser pulse. 6 

In the previous paper we restricted our study to the non-relativistic regime but 
here we consider weakly relativistic effects that occur during gas ionization by a 
strong microwave field. 7 In the present paper we consider the interaction of circularly 
polarized microwave pulse fields with frequency about w ~2x I0 10 s _1 — 2 x lO 11 ^ 1 
with a neutral gas taking into weakly relativistic effects. This pulsed radiation source 
is capable of generating the radiation with an intensity of about l0 8 Wcm~ 2 , whose 
electric field being Eq < lO^Vcm^ 1 is much weaker than the atomic field E a . In 
weakly relativistic case <C c, for example in the weak ionized gas at low plasma 
density, the electron average drift velocity v zav arises in the wave propagation direc- 
tion. At uj ps 2 x I0 10 s _1 and E ps lO^Vcm" 1 we get v E ~ 0.9 x 10 9 cms _1 and 
Vzav = 2.6 x 10 7 cras _1 . If plasma density could get its critical value the drift current 
in plasma ought be j z ps 0.5Acm~ 2 . 8 

In this case, also we study the electron distribution function (EDF) and the 
stability of the discharge plasma in the aforementioned frequency range at weakly 
relativistic electron oscillation energy e e . 

This work is organized in four section. In sec. II. we will obtain produced EDF 
generated by the interaction of the circularly-polarized MW field with a neutral gas. 
In sec. III. we will obtain the dielectric tensor element for this produced plasma. In 
sec. IV. we will study the stability of produced plasma and find the wiebel instability 
growth rate of wiebel instability caused by the anisotropic property of the electron 
distribution function. Finally, a summary and conclusion is presented. 

II. ELECTRON DISTRIBUTION FUNCTION 
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Under condition (1) the thermal velocity of the electrons in a discharge plasma can 
be neglected in comparison to the electron oscillation velocity in the MW radiation 
field. Since the collision frequency is much smaller than the MW field frequency, we 
can ignore the collisional stochastization of the forced electron oscillation as well. 
Furthermore, if plasma density n e (t) produced by the field during gas breakdown, is 
less than the critical density (that is lUq > oo 2 Le = Anne 2 /m) we can neglect the effect 
of the polarization field. Moreover, the plasma density is assumed to be less than the 
neutral gas density n so that the latter can be considered constant. We also suppose 
that the field was adiabatically switched on in the infinite past. Furthermore, we 
can assume that the MW radiation electric field amplitude E is constant during a 
single field period. Therefore, the kinetic equation for the plasma electrons produced 
in the gas breakdown by a strong pulsed field can be written as follows 

-wr+ v.— + e E + -(vx B ) .— = n u ioniz 5( p), (2) 
ot or c op 

where E (£) and B (£) are electric and magnetic fields of a wave propagating along 
the Z-axis, respectively; £ = uot — k r = u (t — z / c) , and u>i on i Z is the ionization prob- 
ability of the gas atoms. Here <5( p) is the delta function of the electron momentum. 
In the weakly-relativistic limit, we can assuming that the electric field depends only 
on time. In the case of MW breakdown, electron-impact ionization is governed by 
the ionization probability uOi on i Z . For a MW discharge Eq. (2) for /o( P, t) is a homo- 
geneous integro-differential equation whose positive eigenvalue 7( E ) determines 
the avalanche ionization constant. However, we can neglect the right-hand side of 
Eq. (2) in the first approximation and calculate the EDF directly by solving the 
Vlasov equation under the following condition, 

loq > 7( E ) (3) 
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where the avalanche ionization constant r y( E ) is to be determined. 8 In this case, 
the condition (3) depends strongly on the neutral gas density and is well satisfied 
at gas pressures of po ~ 10 — 100 Torr. In this approximation, to calculate the 
electron energy distribution function, we assume the field components to be circularly 
polarized 



E x = E sinc^ot, 



E y = E cosc^o^, 



E z = 0, 



(4) 



where, u is the frequency of MW field. Moreover, electric field amplitude E 
describes the slowly varying (over the field period) microwave pulsed envelope, con- 
sider the condition (3) and solving Eq. (2) by characteristic method the equation of 

motion for electrons can be obtained as follows 

dv 7 



ni- 



di 



eE x , 



dv y 

to— = eE y , 



dv z e 

TO— = -{V X E X + VyEy), 



(5) 



dt c 

from this we find solution of the vlasov kinetic equation-characteristics 

v x = -v E {cosLp - cosy? ), 

v y = v E (sin ip- siny? ), 



(6) 



v z = v( 1 ~ ^(V 9 - <Po)), 
where, (p = Ld t is the MW electric field phase and (p = Ld t is the MW electric field 
phase when the electron originates with zero momentum at time to, v E = eEo/muJo 
is the electron oscillatory velocity in an alternating electric field. For the solution of 
the form fo(v,t) = n e (t)fo(v) we obtain 7 
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fo(v x ,v y ,v z ) = S(v x + v E (cosip- cosip ))5(vy - v E (smip- siny9 )) (7) 

x5 (v z - ^[1 - cos(y? - (p )]j ■ 

The function fo(v x ,v y ,v z ) satisfies the normalization condition J d v / ( v) = 1. 
One can show the product of the two last (^-functions in the distribution function (7) 
is independent of ipo. For this reason, by introducing the following notation 

V x = —v y + ve sin uot , V y — v x + v E cos u>ot , 

we obtain 

S(v x + v E (coscu t - cosip ))5(v y - v E (smu t - siny? )) = 77^— S(V± - v E ) , 

ZTTVe 

where 

v ± = v x + v y = v ±+ v e + 2v e(v x cosu t - v y sin u t), v l = v l + v l- 

Therefore, Eq. (7) is reduced to 

fo(V ± , v z ) = ^HV± - v E )S [v z - V \\ - cos(v9 - y? )]) • (8) 

We can see that, the electron distribution function depends on phase field ip and 
electron origination phase <p . Therefore, we must average of electron distribution 
function over the phase field period and origination phase. The projection of the 
phase portrait of the electrons onto the (v x ,v y ) plane is as follows. The electron 
trajectories uniformly cover a circle of radius v E , whose center precesses about the 
origin of the coordinates and describes a circumference of the same radius at a 
rate equal to the microwave frequency. Consequently, the averaging procedure can 
be performed separately for the longitudinal and transverse velocity components. 
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Averaging over the transverse component we obtain 



_ S[v z -^f[l- cos(y? - <p )] ) 
/oK, v z (<p -<p )) = . = = '- (9) 



III. DIELECTRIC TENSOR 

In the previous section we found the distribution function (9) for electrons in 
a discharge plasma is highly anisotropic with respect to the direction of the MW 
radiation field, which, first of all, should result in the onset of the well-known wiebel 
instability. 7 In order to convince ourselves that this conclusion is valid and to find 
the instability growth rate, we turn to the adiabatic approximation, assuming that 
the instability grows faster than the plasma density. 5 ' 9 In this approximation, we can 
use the following dispersion relation for small perturbation 

| k 2 8ij - kikj - -jEij(uj, k) |=0. (10) 

Here Sijiu, k), the dielectric permitivity tensor of the plasma, is obtained by lineariz- 
ing Vlasov equation for the electrons 

9 ~t+ v.§^ + e(E + i( vx Bo)).|£ = 0, (11) 
at a r c op 

and considering self-consistency effect and Maxwell equations and cold ions approx- 
imations as follows 

£i >, k) = 1 + 56^ + 56% = i- ( ^5 ij + ( ^ Jdv 



ViT, VViVj 

ovj uj — k. v 



• (12) 



Here, the second term, i.e., Se^ shows the ion contribution of the plasma dielectric 
permitivity tensor and the third term, i.e., Sefj is related to the electron contribu- 
tion of the plasma dielectric permitivity tensor. The latter kinetic equation for the 



perturbed EDF is valid in the adiabatic approximation when short-wavelength limit 
is justified, i.e., oo Pi <C kv < LO Pe . To solve Eq. (12) by using EDF (9) we introduce 
the following system of coordinate in which k lies only in the XZ plan 



k = (k ± ,0,k z ) , 



Under this condition the non zero contributions of the electrons in the dielectric 
tensor (12) are 



LO 



Pe 



LO 



h 2 

1_ ^ 



2k 2 



2\1 



+ 



(1- 



2(1 



f3 2 UJ 2 ) 



Se 



LO 



.33 



Pe 



2l0 2 



2ik 1 _v%(\ - cos0)" 



ik±VES 







2k z v^(l — cos (f))s k 2 z v%(l — cos0) 2 s 



(13) 



2 „,2 



a; 2 /? 2 



5ei 3 = 



Pe 



2lo 2 



k±v%(l - cos0)g 

A; 2 



a; 2 /? 2 



a; 2 /? 2 



4 (1 



<5ei2 = Se 2 i = 0, 
where 

2 
1 


2 
1 
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LO' 
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for 

for 
for 
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_ 4k 2 v% 

1 w 2 f3 2 

Reu ^ 
Reu = 

ImLO = 
imcj ^ 
Reu = 



' 2lo 2 c 



> when 



1 - cos0) + 5e 13 , 



4A- 2 v 2 



(14) 
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and (3 = (1 - - cos 0) , = ip - <p also, i = 



IV. STABILITY OF THE PRODUCED PLASMA 

As known the process of the gas breakdown by the MW fields is unstable with 
respect to the excitation of longitudinal electric fields and transverse magnetic fields. 
The former is caused due to the positive derivative of the EDF and the latter, treated 
below, is related to the anisotropy of the EDF. 9 In this step we study the weibel in- 
stability, by analyzing electron perturbations propagate across the microwave field. 
Therefore, by averaging dielectric tensor (13) over the <fi between [0, 2n] and sub- 
stituting into the dispersion equation (10) by assuming 2kv 2 E /ujc < 1 yields the 
dispersion relation for electron perturbations propagating across the MW radiation 
field (k ± = 0) 

3 u 2 Pe k 2 v 2 E u 2u 2 Pe k 3 v E 
+ (k 2 c 2 + c4J + c(k 2 c 2 + u 2 Pe ) 1 J 

With assumption k 2 v E <C oj Pe k 2 c 2 and Solving Eq. (15) we obtain the following 
expression, which characterize the growth rate of the weibel instability as follow 

«,= , +i kVEUJPe (l + ^| (16) 



^k 2 C 2 + U 2 Pe ^k 2 C 2 + UJ 2 Pe \ 2w|e 



V. CONCLUSION 

Using a simplest model, we are calculated the weakly-relativistic EDF for plasma 
produced by the interaction of an intense microwave pulse field with a neutral gas. 
The resulting EDF, which can be drive analytically for circularly polarized MW field, 
is highly anisotropic, which indicates that microwave driven plasmas are subject 
to weibel instability. In the non-relativistic produced plasma 7 the electrons only 
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oscillate on the MW field direction but on the weakly-relativistic produced plasma 
the electrons have velocity perpendicular the radiation field. In this case the electrons 
originating with zero energy are entrained by MW field in a certain phase and are 
accelerated to weakly-relativistic velocities, so that the accelerated electrons move 
essentially in phase with the MW field. Then, this process might be repeated for the 
electrons originating over the next time interval, and so on. Analyzing the dispersion 
equation and obtaining the growth rate of the instability produced we found that 
the growth rate modified by an exceed positive term with respect to non-relativistic 
case. Also the instability frequency has a very small real part that can be represent 
the entrain of electrons in MW fields. 
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